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Context. The remarkable stability of extragalactic jets is surprising, given the reasonable possibility of the growth 
of instabilities. In addition, much work in the literature has invoked this possibility in order to explain observed jet 
structures and obtain information from these structures. For example, it was recently shown that the observed helical 
structures in the jet in S5 0836+710 could be associated with helical pressure waves generated by Kelvin-Helmholtz 
instability. 

Aims. Our aim is to resolve the arc-second structure of the jet in the quasar S5 0836+710 and confirm the lack of 
a hot-spot (reverse jet-shock) found by present observing arrays, as this lack implies a loss of jet collimation before 
interaction with the intergalactic medium. 

Methods. In this work, we use an observation performed in 2008 using EVN and MERLIN. The combined data reduction 
has provided a complete image of the object at arc-second scales. 

Results. The lack of a hot-spot in the arc-second radio structure is taken as evidence that the jet losses its collimation 
between the VLBI region and the region of interaction with the ambient medium. 

Conclusions. This result, together with the previous identification of the helical structures in the jet with helical pressure 
waves that grow in amplitude with distance, allow us to conclude that the jet is probably disrupted by the growth of 
Kelvin-Helmholtz instability. This observational evidence confirms that the physical parameters of jets can be extracted 
using the assumption that instability is present in jets and can be the reason for many observed structures. Interestingly, 
the observed jet is classified as a FRII object in terms of its luminosity, but its large-scale morphology does not 
correspond to this classification. The implications of this fact are discussed. 
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1. Introduction 

Jets in Active Galactic Nuclei (AGN) involve some of the 
most energetic processes in the Universe. The way they 
propagate and interact with the ambient medium and, in 
particular, their stability properties have bee n thoroughly 
studied (see Hardee 200fj 120111: IP erucholl201ll . for reviews). 
These jets are subject to the growth of different instabilities 
and many of the structures observed like knots, bendings, 
helices, have been interpreted as a result of this physical 
process. Although there could be other ways to produce 
these structures, such as interactions with clumps of dense 
gas or precession of the central engine, these in turn can also 
give rise to the growth of the instabi lities via coupling to 
any of the unstable mode s (see, e.g. iPerucho et al.l [2005. 
2006 : IMizuno et al l [20071: iMcKinney fc Blandfordl [2 009: 
Mizuno et all 120091: iMignone et al.l 120101: IPerucho et al.l 
2010; IMizuno et al.1 l2011l for numerical studies of cou- 
pling to Kelvin-Helmholtz (KH) and current-driven (CD) 
instabilities in relativistic flows). Our interest in under- 
standing how this process works in jets lies in the pos- 



sibility of obtaining the physical parameters of the flow 
(velocity, gas density and sou nd speed), via modeling 
the observed structures (e.g.. lLobanoy fc Zensud 120011 : 
lHardee et all 120051: lHardee fc Eilekl 120 111) . The Very Long 
Baseline Interferometry (VLBI) technique has provided 
high-resolution observations that produce highly detailed 
images of the structure of AGN jets, and that e ven re- 
solves them transversally (jLobanov fc Zensusll2001ri . Aside 
from those studies in which it is assumed that the observed 
structures are generated by instabilities, no independent 
proof that this is ind eed the case had been provided until 
IPerucho et al.l ()2012f ). where the authors showed that the 
ridge line of the jet in S5 0836+710 behaves as expected if 
it is interpreted as a pressure wave generated by KH insta- 
bility. They also showed, within this framework, that the 
amplitude of the wave grows with distance from the core, 
thus fulfilling the expected characteristics of an instability. 

The luminous quasar S5 0836+710 at a redshift z = 2.16 
hosts a powerful radi o jet extending up to kiloparsec scales 
([Hummel et al.lll992l ). At this redshift, 1 mas ~ 8.4 pc (see 
MOJAVE database and IWrightl l2006t Fl. VLBI monitoring 
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of the source (jOtterbein et al.lli998l ) has yielded estimates 
of the bulk Lorentz factor 7j = 12 and the viewing an- 
gle Oj = 3° of the flow at milliarcsecond scales. The pres- 
ence of an instability developing in the jet is suggested by 
the ki nk structures observed on these scales w i th gr ound 
VLBI (jKrichbaum et al.lll990D . ILobanov et all (|l998j) ob- 
served the source at 5 GHz with VSOF0 and also reported 
oscillations of the ridge-line. Identifying these structures 
with KH modes, they were able to derive an estimate of the 
basic parameters of the flow. High dynamic range VSOP 
and VLBA (Very Long Baseline Array of National Radio 
Astronomy Observatory, USA) observations of 0836+710 at 
1.6 GHz indicated the presenc e of an oscillation at a wave- 
length as long as ~ 100 mas ()Lobanov et al.l [20061 ). which 
ca nnot be readily recon ci led with the jet parameters given 
in ILobanov et all (jl998H . iPerucho fc LobaiKwl (|2007l 120111 ) 
have shown that the presence of a shear layer allows fit- 
ting all the observed oscillation wavelengths within a single 
set of parameters, assuming that they are produced by KH 
instability growing in a cylindrical outflow. In this picture, 
the longest wavelength corresponds to a surface mode grow- 
ing in the outer layers, whereas the shorter wavelengths are 
identified with body mod e s deve loping in the inner radii of 
the jet. In iPerucho et a l. (2012), the authors have verified 
the relation between the helical structure and waves with 
growing amplitude. 

In this paper we present a new observation of the jet 
in this source using EVN+MERLIN, fol lowing the exper- 
iment presented in ([Hummel et al.l fl992f ). We confirm the 
observation of disrupted jet structure at arc-second scales 
and, using our present knowledge of the source, interpret 
this observation. 

2. Observations 

The observations reported here were performed on 1 
March 2008, and lasted 9.5 hours. We used the Multi- 
Elcmcnt Radio-Linked Interferometer Network (MERLIN, 
United Kingdom)0 and, simultaneously, a subset of the 
European VLBI Network (EVN)0. On the one hand, there 
was a total of 9 participating antennas at the EVN: 
Effelsberg (100m diameter), Cambridge (32m), Lovell 
(76 m), Onsala (25 m), Medicina (32 m), Noto (32 m), 
Torun (32 m), Urumqi (25 m), and Shanghai (25 m); on 
the other hand, there were 7 participating antennas at 
MERLIN (some of them simultaneously used as EVN sta- 
tions): Defford (25 m), Cambridge (32 m), Knockin (25 m), 
Darnhall (25 m), Markll (25 x 38 m), Lovell (76 m), and 
Tabley (25 m); 

The observations were performed at 1.6 GHz. The 
recording rate per antenna at the EVN was 512 Mbps (full 
polarization, 8 frequency sub-bands of 8 MHz each, 2-bit 
sampling). The overall synthesized bandwidth at the EVN 

resul ts, and Ned Wrights Cosmology Calculator 
http:/ /www. astro. ucla.edu/~wright/CosmoCalc. html 

2 VLBI Space Observatory Program, a Japanese-led space 
VLBI mission operated in 1997-2007 by the Institute 
of Space and As tronautical Science, Sagamihara, Japan 
http:/ /www.vsop. isas.jaxa.jp/top. html 

6 MERLIN is a National Facility operated by the University 
of Manchester at Jodrell Bank Observatory on behalf of STFC 

4 The European VLBI Network is a joint facility of European, 
Chinese, South African, and other radio astronomy institutes 
funded by their national research councils 
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Fig. 2. Same as Fig. [TJ where a straight line following the 
VLBI jet direction has been drawn to help the eye. 



was 64 MHz, but the total bandwidth of MERLIN was lim- 
ited to 16 MHz. The EVN data were correlated at the Joint 
Institute for VLBI in Europe (JIVE, the Netherlands). 

We used sources 4C39.25 and 3C345 as fringe find- 
ers, and B0859+681 as the amplitude/phase calibrator (3.9 
degrees of separation). The initial amplitude and phase 
calibration was performed with the Astronomical Image 
Processing system (AIPS, NRAC0) using standard proce- 
dures. However, the phases of the different sub-bands in 
the E VN observations were al igned as described in Sect. 
4.1 of lMarti'-Vidal et al.l (|2012ft . due to non-constant sub- 
band phase offsets observed at Effelsberg and Urumqi. 

Once the MERLIN and EVN observations were inde- 
pendently calibrated, the data were exported to the pro- 
gram CASA (Common Astronomy Software Applications, 
NRAO), version 3.3.0, for further calibration, reduction, 
and imaging. We used the ms tool of CASA to scale the 
amplitudes of the EVN visibilities to those of MERLIN at 
similar baseline lengths and imaged the combined dataset 
with the task clean, using natural weighting of the visi- 
bilities and a multi-scale deconvolution algorithm. We used 
beam deconvolution scales of 7 mas (i.e., the size of the EVN 
synthesized beam), 170mas (i.e., the size of the MERLIN 
beam), 50 mas, and 100 mas (i.e., beams of intermediate 
size). We also used the ms tool to apply a phase gradient 
(i.e., a shift in the image plane) of about 1.5 mas (south- 
west direction) to the EVN visibilities, to minimize the im- 
age residuals after the deconvolution (this is probably due 
to effects of source structure, which may introduce small 
offsets in the source intensity peaks at the very different 
resolutions achieved with MERLIN and the EVN). Our fi- 
nal image (with a convolving beam of 84x68 mas; position 
angle of 57 degrees) is shown in Fig. [T] 

The radio map presented in Fig. [TJ shows an irregular 
structure covering a region between 1 and 2 arc-seconds at 
the South- West from the radio core. A large portion of this 
region falls to the West of a straight line extension of the 
inner jet (see Fig. [2]). 



° The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 
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Fig. 1. EVN+MERLIN image of S5 0836+710 at 1.6 GHz. The contours are at 0.12, 0.15, 0.25, 0.33, 0.65, 1.2, 2.5, 5, 
10, 20, 40, 80, and 99% of the image peak intensity (2.48 Jy/bcam). The convolving beam is shown at the bottom-left 
corner. 



3. Discussion 

The large-scale structure is observed at distances from the 
radio core larger than 1 arc-second, after the jet emission 
disappears at 0.2 arc-seconds. This implies an increase in 
the particle energy at the emission site. In typical FRII 
jets this happens where the jet interacts with the ambient 
medium, typically the IGM, and a bright reverse jet-shock 
(hot-spot) is observed in many of these sources. At the hot- 
spot, the collimated jet flow is decelerated and its particles 
are deflected into the radio-lobes. The bow-like shape of 
the radio-lobes resembles the expected forward bow-shock 
in the ambient medium that is difficult to observe. 

In the case of S5 0836+710, the images of the jet at 
the largest scales do not show any signs of enhanced emis- 
sion that could be identified as the hot-spot or any clear 
lobe-like structure. On the contrary, the radio emission 
is irregular. The lack of a reverse jet-shock could imply 
the loss of jet collimat ion, as observed in FRI jets (e.g., 
iPerucho fc MartJ [20071) . After the loss of collimation, the 
jet is transfor med into a subrelativi stic or mildly relativis- 
tic broad flow (jPerucho et al.l 12005! ) that continues propa- 
gating downstream and eventually interacts with the ambi- 
ent medium. Possibly, the observed radio emission is then 



related to the acceleration of particles in the interaction re- 
gion, which, depending on the flow velocity, could involve 
weak shocks. 

In the case of FRI jets, the deceleration and loss of 
collima tion has been related to mass-loading from stellar 
winds flLaing fc B ridle 2002]), strong recollimation shocks 
(IPerucho fc Martd 12007^ 7" growth of nonlinear instabili- 
ties ( Rossi et al.l 120081) or inho mogeneities in the ambi- 
ent medium (jMeliani et alJl2008f) . Observations seem to fa- 
vor the first two possibilities in the region where FRI jets 
are efficiently decelerated, i.e., ~ lkpc. In the case of S5 
0836+710, the VLBI jet is observed to keep collimation over 
more than 100 milli-arc-seconds, which at a viewing angle 
of 3° and z = 2.16 implies ~ 10 — 20kpc. At this distance, 
it seems difficult that mass-loading from stellar winds, even 
if accumulated along the path of the jet through the host 
galaxy, has anything to do with the loss of collimation. In 
addition, there are no signs of strong shocks along the VLBI 
jet. 

Different authors have reported kinks in the jet at dif- 
ferent scales ( s ee the Introduction of this paper) and, in 
IPerucho et al. (|2012l ) the authors identify wave-like pat- 
terns in the observed helical structures. In addition, the am- 
plitude of the helical motion increases with distance. These 
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features are expected when an instability grows in a jet. 
Within this context, a small perturbation to the VLBI jet, 
namely a difference in pressure between two opposite sides 
of the jet, could couple to an unstable Kelvin-Helmholtz 
(KH, if the jet is matter dominated) or current driven (if 
magnetically dominated) mode that grows with distance to 
nonlinear amplitudes, which can lead to jet disruption and 
mixing with the ambient medium. There is a widely held 
belief that jets are dominated by particles beyond the most 
compact region, due to processes like mass-loading from 
stellar winds or clouds, or via the conversion of magnetic 
energy into kinetic energy of the particles in the process 
of jet acceleration (jKomissarov et aD 120091 ) . Thus, the KH 
scenario is favored. 

Linear st ability analysis using expe c ted pa rameters for 
the source (|Perucho fc Lobanovl I2007I . 120111) shows that 
typical minimum growth lengths of disrupting, unstable 
modes such as the surface mode and the first body mode 
are Xi ~ 10 2 i?j. This distance provides the length after 
which the unstable mode amplitude is m ultiplied in an 
expon ential factor. The 1.6 GHz images in iPerucho et al.l 
( 2012]) show that the VLBI jet is observed up to ~ 140 mas 
projected, which corresponds to ~ 22.5 kpc at a viewing 
angle of 3°. The observed growth in amplitude from the 
ridge line of the 1.6 GHz jet changes from 1 — 2 mas at 
z < 50 mas projected (~ 8 kpc deprojected) to 4 mas at 
z ~ 100 — 140 mas projected (~ 16 — 22.5 kpc deprojected). 
Taking into account that the amplitude of the instability 
evolves as A(z) = A exp(z / A,): 



A{z 2 ) = A( Zl ) 



exp(z 2 /Xi) 
exp(zi/Xi)' 



(1) 



For Z\ = 8 kpc, and z 2 = 24 kpc and using a typical 
jet radius of Rj ~ 20 mas ~ 170 p c at 1.6 GHz as ob- 
tained in IPerucho fc Lobanovl (|2007h . A(z 2 ) = 4.7 mas if 
A(z\) — 2 mas. Taking into account that the jet and am- 
bient properties must change in such long distances and 
that this has a direct influence on the growth lengths of 
the unstable modes, the estimate is consistent with the ob- 
servations. This consistency shows that the growth of the 
instability obtained from the linear analysis can explain in 
a self-consistent way the observ ed growth in the amp litude 
of the jet ridge-line at 1.6 GHz (IPerucho et alJl2012j) . 

The growth of the instabilities is possible due to 
the continuous dissipation of kinetic energy from the jet 
flow so mild deceleration is e xpected to occur along the 
observed VLBI jet ( see, e.g. IPerucho et ail 120051 120101: 
lHardee fc Eilekl 201 lh . Deceleration and jet expansion 
would thus be the cause of the gap in emission between 
0.2 and 1 arc-second, within our model. 

The obse r vation al evidence obtained so far 
(jPerucho et al.l I2012L and here) points clearly towards 
this physical process taking place in the jet of S5 
0836+710. Several helical modes, including the more 
disruptive low order body modes, maybe the surface mode, 
and possibly some elliptical modes could be growin g in the 
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jet (|Lobanov et al.lll998l 120061 IPerucho et aLll2012j T These 
wave modes generate helically twisted pressure waves, with 
a maximum that can be relate d to the peak of em ission 
along the jet (the ridge- line, IPerucho et al.l 12012ft . We 
conclude that this FRII jet is probably disrupted by the 
growth of a helical instability. This conclusion represents 
the first strong observational evidence of the development 



of instabilities in jet flows and validates many studies that 
have been made under the assumpti o n tha t jets develop 
such instabilities (e.g.. iHarded 120001 120031 : lHardee et al.l 
120051 : lHardee fc Eilekll2011h . 

In addition, if this FRII jet is indeed disrupted by 
the growth of instabilities, the immediate question is 
whether there are other similar cases. If there are more 
cases, the morphological separation between FRI and 
FRII jets, depending on jet power, should be reconsid- 
ered accordingly. In particular, this could be a possible 
explanation for the different morphology of the jet and 
counter-jets in HYMORS (hyb r id morphology sources, e . g . , 
iGopal-Krishna fc Wiital 12000: iGawrohski etall [2000) . if 
one of the two jets is perturbed and the other is not. 
This requires asymmetric and irregular surrounding me- 
dia (under the assumption that the jet and counter-jet in 
a source are symmetric), but moving clouds in the Broad 
Line and Narrow Line regions can interact with the jet lat- 
erally, inducing internal shocks or sound waves that can 
couple to unstable modes. The observations actually show 
helical structures in th e FRI-like side of some sources (see 
IGawrohski et al.| [2006) . In the case studied here, it has been 
pointed out that the longest helical wavelength observed 
could be produced by long-term precession of the jet noz- 
zle, with periods of the order of Td r — 10 7 yrs, which is a 
similar value to the one given by lHardee et all (|1994l ) for 
3C 449. Future work should be focused on trying to ver- 
ify that the growing instabilities are KH and to obtain as 
much information as possible from the linear perturbation 
theory. 
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